Anomalous magnetic moment of an electron near a dispersive surface 
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Clianges in tlie magnetic moment of an electron near a dielectric or conducting surface due to 
boundary-dependent radiative corrections are investigated. The electromagnetic field is quantized 
by normal mode expansion for a non-dispersive dielectric and an undamped plasma, but the electron 
is described by the Dirac equation without matter-field quantization. Perturbation theory in the 
Dirac equation leads to a general formula for the magnetic moment shift in terms of integrals over 
products of electromagnetic mode functions. In each of the models investigated contour integration 
techniques over a complex wave vector can be used to derive a general formula featuring just integrals 
over TE and TM reflection coefficients of the surface. Analysis of the magnetic moment shift for 
several classes of materials yields markedly different results from the previously considered simplistic 
'perfect reflector' model, due to the inclusion of physically important features of the electromagnetic 
response of the surface like evanescent field modes and dispersion in the material. For a general 
dispersive dielectric surface, the magnetic moment shift of a nearby electron can exceed the previous 
prediction of the perfect-reflector model by several orders of magnitude. 



I. INTRODUCTION 

The coupling of the quantized electromagnetic and 
electron fields to each other gives rise to radiative cor- 
rections in quantum electrodynamics. One of the quan- 
tities altered by these radiative corrections is the elec- 
tron's magnetic moment, which is particularly interest- 
ing because it can be measured to staggeringly high pre- 
cision [1, 2]. The presence of material boundaries affects 
the fluctuations of the electromagnetic field, thus alters 
radiative corrections, and thereby causes the magnetic 
moment for an electron near a surface to differ from its 
value in free space [3-9]. This could in principle be an 
important effect to consider since precision g factor mea- 
surements of leptons not only provide stringent tests of 
quantum electrodynamics but also potentially open up a 
low-energy route to testing physics beyond the standard 
model. 

Previous literature on this subject has shown that mea- 
surement of the shift is not within the reach of contempo- 
rary experiments. However, all these previous investiga- 
tions have made the crude simplification that the surface 
may be regarded as perfectly reflecting. This simplifies 
calculations, but the perfect-reflector model has obvious 
physical deficiencies: it does not account for electromag- 
netic field modes that are evanescent outside the medium, 
neither does it reproduce the fact that any real medium 
becomes transparent at high frequencies. Thus, to make 
a realistic prediction of the potential measurability of this 
effect, one needs to consider a surface which is imper- 
fectly refiecting and dispersive. 

The idea that the use of a more realistic model of 
the material may significantly affect the shift is not sug- 
gested by the results of other electrodynamic boundary- 
dependent effects. For example, in the calculation of 
the Casmir- Polder force on an atom in front of a surface 
[10], the perfect reflector model is completely adequate 
for estimates and its results are reproduced in the ex- 
pected limiting cases, e.g. taking the refractive index of 



FIG. 1. One-loop reducible and irreducible corrections to the 
bare vertex 



a non-dispersive medium to inflnity [11] or the plasma 
frequency of a plasma surface to inflnity [12]. The situa- 
tion is strikingly different for the spin magnetic moment 
of an electron near a surface where, as we will show, 
different models for the electromagnetic response of the 
surface give drastically different results, not necessarily 
obtainable as limiting cases of each other. Something 
similar has already been observed for the self-energy shift 
of a free electron near a non-dispersive surface [13] and 
has since been investigated also for dispersive dielectric 
and conducting surfaces [14]. The reason for these dis- 
agreements between different models has been found to 
be related to the fact that a free particle admits excita- 
tions of arbitrarily low frequency which are dealt with 
very differently in the various models. However, the low- 
frequency part of the photon spectrum turns out to dom- 
inate the shift, so that one necessarily gets completely 
different results for the self-energy in models with dif- 
ferent low-frequency electromagnetic response. By con- 
trast, for bound atomic electrons the gap to the nearest 
energy level provides a natural low-frequency cutoff for 
electronic excitations, so that for the Casimir-Polder shift 
due to a surface any difference in the low-frequency be- 
haviour of the model of the electromagnetic response of 
the surface turns out not to matter to leading order. 

In order to determine the quantum-fleld theoretical 
corrections to the electron's magnetic moment, one usu- 
ally calculates the vertex diagram to the respective order 
of interest (cf. Fig. la to one- loop order). The reducible 
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FIG. 2. Possible n loop insertions to the external legs of 
the vertex diagram. A double line represents a boundary- 
dependent propagator 



one-loop diagrams (cf. Figs. Ib-d) are taken care of 
through mass and charge renormalization. Thus, in free 
space the evaluation of the magnetic moment to one- loop 
order e^ = a is a straightforward calculation and grad- 
uate textbook material (cf. e.g. [15]). However, this 
changes radically if the electromagnetic field is subject to 
reflection by a surface. Then even one-loop calculations 
of quantum electrodynamics get rather complicated [16], 
firstly because of the loss of translation invariance, and 
secondly because of the localization of the electron that 
is required if the calculations are to make physical sense 
and be interpretable. The photon propagator becomes 
boundary-dependent, so that standard mass and charge 
renormalization are no longer applicable. Mass renor- 
malization is affected at all orders, and charge renormal- 
ization from two-loop level up, as illustrated by Fig. 2, 

Most of such technical problems can be avoided by tak- 
ing a different, more appropriate approach. As we seek 
only the correction to the magnetic moment due to the 
presence of the surface, it follows that any boundary- 
independent contributions can be discarded in our cal- 
culation. The bare electron propagator is not affected 
by the presence of boundary, provided the electron is at 
least a few Compton wavelengths away from the surface 
of the medium so that the interaction between them is 
wholly electromagnetic, which is the case we are inter- 
ested in here. Then the one-loop vacuum polarization 
diagram can be discarded when coupled to a static exter- 
nal field, and to one-loop order the charge and its renor- 
malization are the same as in free space. None of the 
three remaining one-loop diagrams contain any particle- 
antiparticle loops. Consequently, for the calculation of 
the boundary- dependent corrections to the magnetic mo- 
ment at one-loop level it suffices to work with a first- 
quantized electron interacting with a second-quantized 
photon field. This allows us to borrow well-tried tech- 
niques from quantum optics. We shall use perturbation 
theory to determine the energy shift — /x • Bq, due to the 
presence of the material surface, of the electron's spin in 
a weak external magnetic field Bq, and then extract the 
magnetic moment as the coefficient fi of the term linear 
in Bq. 

In the next section, we give the derivation of the shift 
in terms of mode functions of the electromagnetic field. 
Then, Sec. Ill A deals with the electromagnetic mode 
functions for a non-dispersive dielectric. Sec. IIIB with 



those for a plasma model, and Sec. Ill C with the case of 
a dispersive dielectric. In Sec. IV we show in detail how 
to evaluate the required integrals for the various differ- 
ent models, and in Sec. V we discuss the results of the 
various models. Finally, in Sec. VI we discuss the exper- 
imental relevance of our results. We use natural units 
with c — I ^ h and eq = 1 = Mo throughout. 



II. PERTURBATION THEORY 

Our starting point is the Dirac equation coupled to an 
electromagnetic field A^ 



^il^{df, + ieA^) + m]?/; = 0. 



(1) 



In order to best exploit the similarities with quantum 
optical problems we shall use the Dirac equation in its 
non-covariant form 



d 
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(2) 



where 7° = /3, 7* = /3a' and A^ = {$,— A). In or- 
der to derive the energy shift of a stationary electron in 
a weak and static magnetic field and close to a reficct- 
ing surface, we aim at a non-relativistic approximation. 
A routine method of achieving this is to approximate 
the Dirac equation by means of a Foldy-Wouthuysen 
transformation, as has been done for working out the 
magnetic-moment correction of the electron near a per- 
fect reflector [9]. However, this approach requires great 
care in all its steps since several successive orders in the 
non-relativistic expansion turn out to contribute to the 
shift. Additionally, in the only comparable previous lit- 
erature [9] a second unitary transformation is applied to 
the transformed Hamiltonian (an approach that the au- 
thors term the 'Paris Method'). This second transforma- 
tion only works for the perfect reflector, meaning that 
a straightforward extension of Ref. [9]'s work to other 
surfaces is not possible. Hence we adopt a more gen- 
erally applicable approach. We work directly with the 
Dirac equation and apply perturbation theory using the 
solutions of the Dirac equation in a static magnetic field 
introduced in Ref. [17] as unperturbed 4-spinor states. 

We subject the electron to a constant classical mag- 
netic field Bq acting initially the along z axis; Bq = B^z^ 
as shown in Fig. 3. Generalization of our results to in- 
clude the case of magnetic fields parallel to the surface 
will be provided later. A suitable classical vector poten- 
tial that generates such a field Bq is Aq = —5(1" x Bq). 
Added to this is the quantized photon field Ag, so that 
the total vector field entering the Dirac equation (2) is 



Ao 



^Q 



(3) 



The Dirac equation for a particle in a constant mag- 
netic field can be solved exactly, whence the unperturbed 
Hamiltonian is 



^0 



Pm 



with 



eAo 



(4) 
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FIG. 3. Basic setup 



We show in Appendix A how the eigenstates of the 
Schrodinger Hamiltonian for a particle in a constant mag- 
netic field can be obtained and how one can easily gen- 
erate those of the corresponding Dirac Hamiltonian from 
them. The Hamiltonian can be expressed in terms of 
the creation and annihilation operators for right circular 
quanta that serve to move states from one Landau level 
to the next and obey the bosonic commutation relations 



[Sfl,s^] = i, [bRM = Q = [bWb^R] 



(5) 



The perturbation (Dirac) Hamiltonian that describes 
the interaction of the electron with the quantized field is 



Hu 



-e7o7 • Aq = —ea ■ A, 



Q 



(6) 



We note that the scalar component $ of the quantized 
field shifts all states uniformly, so has no impact on the 
magnetic moment. The energy shift is given by second 
order perturbation theory as: 



A^ 



T 



d-^k^ 



E - E' 



(7) 



where Ik. a indicates a one-photon state with wave vec- 
tor k and polarization A, and ^e represents the state of 
the electron coupled to the classical field Aq. The Dirac 
eigenstates ^e for each Landau level may be written in 
terms of the corresponding Schrodinger eigenstates {i'), 
as shown in Appendix A. 

The quantized electromagnetic field is written in terms 
of mode functions fkA 



Aq = ^ / d^k (fkAflkAC" 






e*"*) (8) 



:;t 



where dkA a-nd d|.^ arc the photon annihilation and cre- 
ation operators for each wave vector and polarization. 
Any normalization constants that may appear for the 
modes in different classes of models have been absorbed 
into the functions fkA in order to preserve the generality 
of the expressions to be derived. 

Substituting the quantized vector field Aq from 
Eq. (8) and the Dirac eigenstate ^e as derived in Ap- 
pendix A into Eq. (7), we have for the energy shift of an 
electron in Landau level v and spin state s: 



AE 



=^EE 



d-^k 



I {i^, s\ {Ho + E,)a ■ fkA(go + E,,) \v' , s') p 
A{E^ - E^. - uj)E^[E^ + m)\E^.\{\E^, \ + m) 



(9) 



where w is the photon frequency. As we are aiming for 
the energy level shift of a localized stationary or slow- 
moving electron, we perform a non-relativistic expansion 
which implies that the electron's rest mass is much big- 
ger than all other relevant energy scales. In this sense 
we expand, loosely speaking, in powers of 1/m, and we 
shall keep contributions to the magnetic moment up to 
order l/rrv'. For the surface-dependent corrections to 
the magnetic moment that we are interested in, this is 
in fact an expansion in l/{mz), which is the ratio of 
the electron's Compton wavelength to the distance z of 
the electron from the surface, and thus very small in- 
deed for all even remotely realizable values of the dis- 
tance z. The above expression contains terms that are of 
order 1/m? and higher; however, it turns out that none of 
the 1/m? terms contribute to the magnetic moment, all 
the contributions to which are of order I/itv" and higher. 
We note that fourth-order perturbation theory can con- 
tribute only terms of order l/m^ or higher, and therefore 



r 



does not need to be considered here. 

Since the spin magnetic moment is obtained from the 
coefficient of the terms of A£' that are linear in (TzEq, one 
needs to carefully enumerate all the possible effects which 
may generate a dependence on Bq. With this in mind, we 
find that we must go beyond the dipole approximation 
for the field Aq, because the multipole expansion 



AQ(r) = AQ(ro) + [(r - Tq) • V]AQ(ro) 



(10) 



ro, 



contains powers of the displacement operator r - 

— 1/2 

whose matrix elements are proportional to Bq (see 
Eq. (A15) in Appendix A). For almost all contributions 
to the energy shift we shall find that the dipole approxi- 
mation works and we need to take along only AQ(ro), but 
there is going to be just one exception for which we shall 
have to take along the next term in the above expansion. 
This may seem unusual from a view point of quantum 
optics, but this is ultimately due to the entirely classical 
effect of the magnetic field causing particle trajectories 



to bend. Wc shall see that to order l/vn? there is only 
one term for which the curvature of the trajectory mat- 
ters and the variation of the quantum field Ag along the 
curved trajectory causes transitions between one Landau 
level and the next. 



Despite us being interested in the energy shifts only for 
particles, not anti-particles, we need to sum in Eqs. (7) 
and (9) over all intermediate states, both particle and 
anti-particle ones. It is convenient to parse the prob- 
lem and consider the arising matrix elements separately 
for particle-particle and particle-antiparticle transitions. 
Using the explicit matrix form of the unperturbed Hamil- 
tonian Hq from Eq. (4) and from now on suppressing the 
subscript kA on the mode functions fk^ that signals its 
dependence on photon wave vector and polarization, we 
have: 



(/fo + E,,){cx ■ f)(iJo + E,)=\ ^^'; ^' 



The explicit forms of the matrix elements are: 

Hee = {E^' + m){(T ■ f)(cr • tt) + {E^ + m){a- ■ 7r)(cr • f) 
i/ee = {o- ■ 7r)(cr . f)(cr • tt) + {E^, - m){E^ + m){(T ■ f) 
Hse = {E^ - m){(T ■ tt){(t ■ f) + {E^, - m)(cr • f)(cr • tt) , 

where the momentum of the particle moving in the ex- 
ternal magnetic field has been denoted by 



eAo 



(12) 



We list Hee only for completeness; it is irrelevant for the 
energy shift of initial states that are particle rather than 
anti-particle states. 

In order to find magnetic moment contributions up to 
order l/m^, we expand each term in Eq. (9) for large m 
and extract all terms proportional to ctzEq, while taking 
care to include, where necessary, the effects of the multi- 
pole expansion of the field, Eq. (10). When manipulating 
the above matrix elements one needs to remember that 
TT contains a differential operator that does not commute 
with any quantity that is a function of position; its com- 
mutator with a function gir) is [7ri,(7(r)] = —i\/ig{r). 
With that in mind wc use the algebra of the a matrices, 



(5,o + ' 



^ijk 



ak , and decompose iJge along the lines 



(11) of Eq. (A14) to rewrite it as 



J 



H, 



il3o{E, + S.'+2m)/_ - {E,, - E,)Po{a- x f)_ 
+ {E, + E,, + 2m)f,p, + {E,, 



- iPo{E, + E,, + 2m)/+ + {E,, - E,)Po{(t x f)+ 



E^){a X f), + {E^, + m)(T • (V x f) 



(13) 



r 



The expression we need to consider is 



^^EE 



d^k 



{v.s\ Hp 



4{E, - E,, - ij)E,{E, 



m)E^f (E^, + m) 

(14) 
Carrying out the large m expansion and discarding terms 
higher than order \/m?, we find that none of the cr x f 
terms contribute to the energy shift to this order, as the 
denominator of Eq. (14) is already yielding terms of or- 
der \/m^. The weak magnetic field expansion i?o — J- of 
the energy difference Ei, — Ey' — uj generates factors of 
cBq / {muj) . Thus the significant terms from the numera- 
tor have to be of order m? . Of those, the terms in f-fl_ 
and f+fX ^^d those involving p^ are spin-independent, so 
do not contribute to the magnetic moment. This leaves 

I 



three terms, whose matrix elements are proportional to 

(z/, s\ /_ \v', s') {v\ s'\ cr • (V X f*) \v, s) + H.C. 
{v, s\ f+ \v\ s') {v', s'\ ^ • (V X f*) \v, s) + H.C. 



( i/, s| cr. (Vxf)|iy',.s') I' + H.C. 



(15) 



One might think that the first two terms cannot con- 
tribute, since the matrix elements on the left cannot 
change the spin s and the matrix elements on the right 
at first sight appear not to be capable of changing v, 
but the latter is true only in the dipole approximation. 
Taking along the next-to-leading term in the multipole 
expansion of cr • (V x f *), we see that they in fact do ad- 
mit transitions between Landau levels. The contribution 
to the energy shift from transitions into the intermediate 
state with v' = v — 1 \s 



Ai?quad(j^' = ^-1) = ^E/ d'\^'-^V^f- {y - l|(r - ro) • V(^ • (V x f*)) \v) + H.C. , 



(16) 



and the contribution from those with v' = v + \ is very similar. Substituting the matrix elements of the position 
operator from Eq. (A15), wc see that the factors of /3o cancel and this contribution is indeed of order B^ in the external 



magnetic field. The sum of the i/' = v + 1 and i'' — ly — I terms can be written as a set of second partial derivatives 
of mode functions. We can simplify the resulting expression if we take into account the orthogonality properties of 
the polarization vectors of the electromagnetic field (which we shall discuss in detail in the next section) and the fact 
that some terms arc odd under the even integration over all photon modes and hence drop out. We find for the total 
contribution from these quadrupole terms 



AE, 



quad 



4m3 



S/^'": 



1 






r ^ J X I r ^ J X r 



d- 



dxdy 



dy"^ 



''Il\ 



dxdy j 



H.C. 



r 



The third line of Eqs. (15) contributes for s ^ s', and in 
this case the dipolc approximation is enough to deliver 
all terms that are relevant for the magnetic moment to 
order 1/m^. This spin-dependent term contributes 



AK, 






^y'rf3kj_(|(vxf)j2+i(vxf)j^ 



(17) 
The sum of quadrupole and spin-dependent contributions 
gives the particle-particle portion of the shift, 



AE, 



Ai?quad + AEt 



spin 



(18) 



We now turn our attention to the matrix element 
between initial particle and intermediate anti-particle 
states, Hee- Disentangling a product of three a matri- 
ces is a slightly lengthier procedure than dealing with 
just two of them in -ffee, but the calculation runs along 
exactly the same lines: one repeatedly applies the algebra 
of the a matrices. In addition to the non-commutation of 
the momentum tt with any function of position, one now 
also has to take note that the vector product tt x tt does 
not vanish but has a z component (tt x 7r)z = — 2i/3Q, 
which comes about because the x and y components of 
TT do not commute, i.e. [7ri:,7rj,] = — 2i/3g. Then straight- 
forward calculation yields 



Fee = - 2/32(S^5y_a_ + bRbRUa+) + blPo {2ipM-^z + <^-fz) + *(V x f)_ - [a x {V x f)]_} 

+ 5fl/3o {2zp.(/+a, + (7+/,) - z(V x f)+ + [^ x (V x f)] + } + 2(5fblbj,{f-a+ + f+a^) + 2filf+a. 
- {a ■ f)7r2 -h ef • Bo + iea ■ (f x Bq) + (cr ■ i){E,,, - m){E„ + m) 

The denominator of Eq. (9) now gives terms of order I/tti^, and non-relativistic expansions of the various energy 
terms supplement this with factors of eBo/rn?. Therefore an obvious contribution to the energy shift of order eBo/mP' 
comes from final term in H^e and whose square contributes with m^ in the numerator. It gives: 



e-'Bn 



AE,,,, ^-^Y^Y. Sk{l + 2u + s + s')\ (,s| a ■ f |,s') |^ 



2m^ 



E/'^'^i/- 



(19) 



The only other possible contributions are those for which the factor of 1/ni^ supplied by the denominator is multiplied 
by a numerator of order P^m? = —eBt)m?/2. Such contributions arise from the cross multiplication of the last term 
in iJee with terms that carry /Jq or eB^. Among those, the terms with b\fbj^ and 6^6^ can be ruled out as they lead 
into intermediate states with v' = u ±2 for which the other factor, the matrix elements of the final term in i^ee; 
vanishes. Therefore the only contribution wc still need to consider comes from a term proportional to 

Y {s\ 2lilv{J_a+ + /+CT_) + 2(ilJ+a- + eBo{(T ■ i){2v + 1) + ef • Bq + tea ■ (f x Bq) \s') {s'\ a ■ f* \s) + H.C. 



As nothing inside the matrix elements depends on s', the 
sum of \s') (s'l gives the identity operator. Multiplying 
out all the a matrices and collecting terms proportional 
to Uz is then straightforward and gives; 



AEee2 



'Baa 



Am^ 



A 



d'kl/,1 



Here we have again taken into account the orthogonality 
properties of the electromagnetic polarization vectors, on 
account of which e.g. J2\I'^'^^ fxfy = 0- Going beyond 
the dipolc approximation is not necessary here and does 
not generate any additional contributions from H^e to 
the accuracy of the energy shift that interests us. Thus 
the energy shift due to intermediate antiparticle states 



is: 



^Eee — Ai^^ee 1 + Ai?eg.2 — 






\ -J 



(20) 
From the total energy shift A_E = AE'ee + Ai^gg, we 
extract foUowing magnetic moment shift, via /S.E = 



A/i_L 






l(Vxf) 



all 
modes 







h 






H.C. 



r 



(21) 



This constitutes the first major result of this paper — an 
expression that delivers the magnetic moment shift for 
an electron interacting with any quantized electromag- 
netic field subjected to boundary conditions. We have 
presented the calculation with an externally applied mag- 
netic field Bo that is perpendicular to the surface of the 
medium in mind, but the orientation of any surface has 
not actually entered the calculation at any stage. Hence 
we can quite easily find the corresponding shift for a mag- 
netic field parallel to the surface by choosing the mode 
functions f (r) such that they are solutions of the wave 
equation with a differently oriented surface of the dielec- 
tric, or equivalently cycle the Cartesian coordinates in 
Eq. (21). One can also apply the formula to the case of 
non-planar surfaces, provided one can convert the mode 
functions to Cartesian coordinates. 

We note that our derivation and the resulting Eq. (21) 
can of course not be used to calculate the anomalous 
magnetic moment in free space. While one could fudge 
an estimate by cutting off the integral over photon fre- 
quencies at w ~ m, which would give the correct order 
of magnitude e'^/m, a correct calculation would require 
the field quantization not just of the photon but also of 
the electron. As explained earlier, if we are interested 
only in the shift of the magnetic moment compared to 
its free-space value and if the interaction with the sur- 
face is wholly electromagnetic so that the electron is not 
subject to boundary conditions due to the presence of the 
surface, the electron propagator is not affected by the in- 
teraction with the surface and we can calculate the shift 
by using a quantum-optical approach and quantize only 
the photon field. 



III. ELECTROMAGNETIC MODE FUNCTIONS 

A. Non-dispersive dielectric 

We wish to determine the magnetic moment shift (21) 
for a range of models for the material and start by consid- 
ering a semi-infinite slab of non-magnetic, non-dispersive 



material that fills the half-space 2; > as shown in Fig. 3. 
The dielectric function is 



e(r) = l + e(z)(?i2_i) 



(22) 



where n^ > 1 is the index of refraction. Maxwell's equa- 
tions in a dielectric but non-magnetic medium with the 
dielectric function e(r) give the following equation for the 
quantized vector field Aq, 



V X (V X Aq) = -e(r) 






(23) 



For a piecewise constant dielectric function, like the one 
in Eq. (22), the solutions of the above differential equa- 
tion are found most efficiently by employing the general- 
ized Coulomb gauge, defined by 



(r)A 



QJ 







(24) 



so that, except directly at the interface z = 0, Eq. (23) 
reduces to the standard wave equation with plane wave 
solutions. These plane wave solutions are then patched 
together at the interface 2; = by imposing the continuity 
conditions that follow directly from Maxwell's equations. 



E||(r)U=o- 
B(r)U.o- 



e(z)i?,(r)U^o+ 
E||(r)U=o+ 
B(r)U=o+ • 



(25) 



In the following we shall write wave vectors that exist on 
the vacuum side as k, and those on the dielectric side as 
k''. A superscript R will denote a reflected wave vector 
with a reversed z component. It will be convenient de- 
compose wave vectors k into components parallel to the 
surface (k||) and perpendicular to it (fcz). Modes will be 
labelled by their direction of incidence, incident momen- 
tum k, and polarization A, and will be broken down into 
incident, reflected and transmitted parts. Our notation is 
such that sgn(fcz) = sgn(fcf) on the real axis, which cor- 
responds to connecting incident and transmitted waves. 
Modes that are incident from inside the dielectric may 
suffer total internal reflection at the interface, and thus 



be evanescent on the vacuum side. This means that there 



is a certain range of values for kf whose corresponding 
kz are pure imaginary. The mode functions are 



fcA.nondisp 



pright 



(2^) 
1 



feA.nondisp (27r)3/2 ^/2^ , 



\ ^ {9(-z)[e""-e,(k) + 



jLiWr- 



I2uj 
1 1 



i?J:e*''"-^eA(k«)]+e(z)ri^e''^--eA 



i?„ik"«r 



i?fe 



eA(k'*^)] + e(-z)ri^e*'''-eA(k)} 



(26) 



r 



where the polarization vectors e\ (k) obey the gauge con- 
dition k • e;^(k) = 0. A convenient choice is a decompo- 
sition into transverse electric (TE) and transverse mag- 
netic (TM) modes through 



eTE(k) = — (fc.y,-fcj;,0) 

fcy 

eTM(k) = ■— f fc^rfc^, kyk^, 



(27) 



tion. They are normalized such that the radiation Hamil- 
tonian is represented as a collection of harmonic oscilla- 
tors, 



ifrad = ^ / t^'r 



A| + (V X Aq)' 



(e(r) 
Y^ J d^k uj UkAflL + 2 ) 



(29) 
(30) 



The reflection and transmission amplitudes in the mode 
functions (26) are given by the Fresnel coefficients 



i?.^ 



/c?- 



TE 



i?XM ~ 



k^ + k'^ 
n^k^ - kf 
ri^kz + fcj? 



rpL 



2A:. 



Tr 



kz + kd 
2nk, 



TM 



n 
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Rx = -Rk 



yd 

rpR _ '^Z rpL 



with z components of the wave vectors in vacuum and in 
the medium connected by the laws of refraction, 



\Jn^{kl 



kl 



(28) 



The modes are the same as in Refs. [13, 16, 18] but with 
slightly different conventions concerning their normaliza- 



The modes form an orthogonal and complete set, as ex- 
plicitly shown in Refs. [18, 19]. That this must be so can 
easily be seen by writing Eq. (23) for the mode functions 
in the form [20] 



1 



Mr) 



zVx 



1 



Mr) 



V X ^^^f7{v)f^x{v) = c.'2vMr)fkA(r) . 

.(31) 

Evidently, this is an eigenvalue problem for a Hermitean 
operator acting on •\/e(r)fkA, which must therefore form 
a complete set of orthogonal functions satisfying 

jd^v 6(r)/*,(r)/k.A'(r) = ^5xy5^^\]^ k) . (32) 

Substituting these modes into Eq. (21), we find that the 
magnetic moment shift can be written in the form: 



A/iJ.^:- 



4m3 



)?=±1,A 



E h'^ 



dfc,g^(fcjl,fc,)[l + |i?tl1 + 



R\2 



dktgim\,kz)\T^\ 



+ d dkzgt{k^^,kz)Rk{e''''^^+e-''''n + - dkfgt{k^^,kz)\T^\ 



nR\2 2ik^z 



(33) 



where the functions g^(fc||,fcz) are 
kl 



5te 



5tm 



/fc3 ' 

2K|| -I- k^ 
2fc3 



5te 



2fc3 



ffTM — 



(34) 
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The critical value of fcf, below which the modes are 
evanescent on the vacuum side, is F = vr?^^ k» . 
The functions ^^(fcii,^^) have a branch cut due to k = 

\ k1 + fcn in their denominators, but are otherwise ana- 
lytic in kz- We place this cut at kz = ztifcii . . . ± ioo. The 



Im(*;,) 



c'( 



-> > Re(fc,) 



-ik] 



Vn^-l 
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FIG. 4. Complex fcz plane for non-dispersive dielectric 



contributions from the various types of modes to Eq. (33) 
are easily identifiable: the integrals over k^ = —oo . . . — T 
and kz ~ . . . oo correspond to right- and left-incident 
travelling modes, respectively, while the integral over 
fcf = —T ... corresponds to evanescent modes, fc^ has 
branch points at k^ = ztiF/n; we place the branch cut 
in between. Using dkf = n'^[kz/kf) dkz, we can manipu- 
late the kz and fcf integrals from the first line of Eq. (33), 
which are all independent of z, into 



dkzgt{ki\,kz 



1 + \R 



L\2 



kz 



IT: 



R\2 



dkzgt{k\\,kz) , 



(35) 



where the equahty follows since kz and fcf are here both 
real. The second line of Eq. (33) contains the z dependent 



terms and can be written as: 

/•oo 

f-iV j n 1 

+{) dkzj^gt{k^^,kz)\T,Ye'"'^' . (36) 

Jo '^z 

We observe that for real kf and pure imaginary kz 

kz 



R^ \kd — ^K ~ R„ 



^^ rriRrjiR:^ I 



(37) 



\k<i = ~K — -n'cr\k<i=K - Ld^'T -CT 

which permits us to combine the integrals (36) into one. 



d / dkzgt{k\\,kz)Ri 
Jc 



L 2ikzZ 



(38) 



with the contour C as shown in Fig. 4. Thus, we have 
for the magnetic moment shift. 



A/ij 



i?=±l,A 



d^ki. 



dkzgl{k\\,kz)Rie 



2ik, 



dkzgt{k\\,kz) 



(39) 



The second term in the brackets can be left out since it 
does not depend on the position z of the particle. It arises 
from the same electromagnetic field fluctuations but in 
vacuum, i.e. without any dielectric medium present, as is 
easy to see by taking the limit n — >■ 1 and hence i?^ — ?> 
in the above equation. Therefore, we subtract it as a 
free-space counterterm and obtain for the renormalized, 
position-dependent magnetic moment shift 



Afij 



— y 



i5=±l,A 



d^kii 



dkzdgt{k\\,kz)R^e 



2ikzZ 



(40) 

The integrand is analytic in the lower half of the complex 
kz plane, so we can freely deform the contour C into that 
labelled as C" in Fig. 4. Carrying out the sum over t? 
and the angular part of the k|| integration, we find the 
following explicit expressions for the magnetic moment 
shifts Afi± and A/i|| in magnetic fields perpendicular and 
parallel to the surface of the dielectric medium: 



A/ij 



A^ll 



327r277l3 



dku 



dk\\ 



dk. 



dkz 



fc3 



2fc3 



(2fc| - kl) Rk^ + (2fc2 + ki) 



k^,]R^ 



TM 



^2ikzZ 



(Sfcf + 2kfj i?^E + (Sfcfi - 2kfj i?^M 



^2ik^z 



(41) 
(42) 



The shift A/i|| for the case of the externally applied mag- 
netic field directed parallel to the surface of the dielectric 
has been evaluated in exactly the same way as A^^ , but 
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by cycling the Cartesian indices in Eq. (21) so as to de- 
scribe a differently oriented surface, as explained below 
Eq. (21). 



The integrals in Eq. (42) can be calculated exactly. 
However, we postpone discussion of their calculation and 
results until after consideration of a plasma surface, in 
order to be able to compare and contrast the two models. 



B. Plasma Surface 

We now consider the simplest model of a dispersive 
medium — an undamped plasma of freely moving charge 
carriers. The dielectric function for a half-space is then 



e(r,a;) 



i-e(.)-| 



(43) 



where ujp = ^ Ne^ jmc is the plasma frequency which 
characterises the material made up from charge carriers 
of effective mass ttIc, charge e, and an average density N . 
As discussed in detail in Ref. [21], Maxwell's equations 
yield the equation to be satisfied by the quantized vector 
field, and thus by the mode functions, as 



Vx [VxfkA(r)]+e(z)c^2f^A(r) 



'fkA(r) . (44) 



Using Coulomb gauge one can simplify this to the wave 
equation, with a shifted frequency inside the material, 
and its solutions are travelling and evanescent waves. 
The continuity of the parallel wave vector k|| means that 
the relation between the z component of the wave vector 
inside the material, fc^, and the one in vacuum, fcz, is 
now given by 

i ■ (45) 



^d2 ^ ^2 



As Wp is real, it is possible for fc^ and fc^ to be simulta- 
neously imaginary, i.e. there exist solutions to the wave 
equation which fall away exponentially on both sides of 
the interface; these are surface plasmons. Eq. (45) also 
tells us that if /cf is real then k^ must also be real, so 
that one cannot have modes originating from within the 
medium that become evanescent on the vacuum side. 
Thus, we expect to reproduce some of the features of 
the perfect reflector model. However, for fc^ < ojp the 
wave vector inside the medium, fc^, is pure imaginary 
and therefore the modes are travelling waves on the vac- 
uum side but evanescent inside the medium. So, there 
are three types of mode: TE, TM and surface plasmon 
(SP). The TE and TM modes are always travelling waves 
in vacuum and can be travelling or evanescent waves in 
the medium. They can be written in almost the same 
form as in Eq. (26) for the non-dispersive dielectric. In 
order to see what is different, we note that the left hand 
side of Eq. (44) is a Hcrmitian operator acting on fkA, 
which implies the orthonormality relation 



y'd^rf^,,,(r)fkA(r) = i-5,,,5(-^)(k' - k) . 



(46) 



The difference between this and the orthonormality rela- 
tion (32) for the non-dispersive case, is a factor of e(r) un- 
der the integral. However, this is compensated by the dif- 
ferent relation between fc^ and fcf which now gives dk'^ = 



{k,Jkf)dk^ and therefore 5{k^ - k'^) = (k^/kf)5{kf - fcf ) 
in Eq. (46). We find that, apart from the different dielec- 
tric function to be used within the reflection and trans- 
mission coefficients, the only difference in the TE and TM 
modes for the plasma as compared to the non-dispersive 
dielectric is an overall factor of n in the right incident 
modes. Thus, the TE and TM modes for the plasma 
are obtained from the non-dispersive dielectric modes in 
Eq. (26) via the replacements 



rlcft 
kAplasma 

fright 
kA, plasma 



cright t 2 



m 



kA.i 



.disp('^ -> <^)) ■ 



The SP mode is derived by considering modes that 
exponentially decay on both sides of the surface. Defining 



ifc. 



kl 



-ikt 



fc2 - e(a;sp)a;2p 



one finds the solutions of Eq. (44) to be of the form 
exp(K2 -I- ik|| • r||) for z < and exp(— k'^z -|- ik|| • r||) 
for z > 0, with the polarization vectors as in Eq. (27). 
The continuity conditions (25) then relate the mode am- 
plitudes on both sides of the interface. For the TE po- 
larization the resulting equations turn out to contradict 
each other whence there is no TE surface plasmon. For 
the TM polarization the combination of these continuity 
conditions yields the relation 



-e{(jJsp)K, 



(47) 



which shows that surface plasmons can occur only for 
frequencies where e(cj) is negative and delivers the surface 
plasmon dispersion relation. 



sp 



^i-^ 



(48) 



As before, we find normalization constants by requiring 
that for each of the TE, TM and SP modes the Hamil- 
tonian reduces to that of a set of harmonic oscillators, 
Eq. (30). For the SP mode the set of oscillators exists 
only for the TM polarization and the sum over modes 
is just a two-dimensional integral over kii , as the disper- 
sion relation, Eq. (48), fixes lo at the surface plasmon 
frequency Wgp , which also fixes k^ . 

The mode function for the surface plasmon is specific 
to the plasma model; it reads [12, 21] 



fk.SP = 



1 



1 



-e(z) (k|| + ^z 



with the norming function 

£4-1 



ik 



e(-z) (k|| Jiz I e 



■ikii -rii -\-Kz 



pih) 



'7^(1+^ 



■rn —K z 



e = e(a;sp) 



(49) 



(50) 
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Substituting the mode functions into Eq. (21), we find 
that, similarly to the case of the non-dispersive dielectric, 

I 



the magnetic moment shift due to the TE and TM modes 
may be written as 



A TE,TM 



— F 



i3=±l,A 



rf^kii 



dk 






where the functions 5^ (fc|| , fcz) are the same as in Eq. (34) 
for the non-dispersive dielectric. 

Again, the z independent terms of this integral make 
up the free-space counterterm to be subtracted. In this 
case, these are all terms except the one proportional to 
i?^. In order to see that these really give the free-space 
contribution, 
as 



we use d/cf = {kz/k'^)dkz to rewrite them 



dkzgt{k\^,kz){l + \RJ: 



L\2\ 



" LJ-n \ Z / 



(52) 



Here kz is always real, as there are no evanescent modes 
on the vacuum side. However, fc^ is imaginary for kz < 
ujp , which means that in the first integral | i?^ P = 1 ■ The 
second integral can be simplified by applying Eq. (35) as 
kz and kf are both real. Therefore, Eq. (52) is just a con- 
stant, independent of the properties of the material and 
the same as one would get in free space. Subtracting it as 
free-space counterterm from A^^^ ' , we arrive at the 
same expression as in Eq. (40), except with the contour 
running straight along the real axis as shown in Fig. 5. 
The integrands (up to the explicit form of e within R^) 
for the TE and TM modes in the plasma model are identi- 
cal to those for the non-dispersive dielectric. Just as with 
the non-dispersive case, we deform this contour into the 
lower half-plane. We note that for the plasma model the 

TM reflection coefficient has a pole at —iK. = \ ^Ip — ^m ■ 



(51) 
I 

We are left with a contribution from the contour inte- 
gral over C" that is identical to Eqs. (41) and (42) (the 
corresponding quantities in the non-dispersive case) , and 
a contribution proportional to the residue of the TM 
reflection coefficient. Now we add the surface plasmon 
part, which is derived by substituting the surface plas- 
mon mode functions, Eq. (49), into the expression (21) 
for A/ij^. This gives the magnetic moment shift due to 
the surface plasmon as 



A^f = 






dk\\k\\ 



2kl 



dkuku 



p(fc||)K2 

3fc| -I- 2k2 
2p(fc||)^2 ' 



(53) 

(54) 



rrrrp 'T'lVT 

Adding these to the results for A/Xj_ ' we find that 
they exactly cancel the contributions from the residue 
part of the contour integration along C", similarly to 
what has been found in the calculation of atomic en- 
ergy level shifts near a plasma surface [12]. Thus for the 
plasma surface the total magnetic moment shift, i.e. the 
sum over TE, TM and surface plasmons contributions, 
is given precisely by Eqs. (41) and (42). In other words, 
the magnetic moment shift for both the non-dispersive 
dielectric surface and the plasma surface can be found 
from the same integral simply by inserting the appropri- 
ate dielectric function into the reflection coefficients. The 
additional modes arising in the plasma model are auto- 
matically taken care of by the more complicated structure 
of the expression in the complex kz plane over which the 
contour is deformed, with the same end result. 
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FIG. 5. Complex kz plane for plasma surface 



C. Dispersive Dielectric 

The next model we would like to consider is a disper- 
sive dielectric. For this we need to move the pole in the 
dielectric function of a plasma away from zero frequency 
to a finite transverse optical phonon resonance lot, which 
corresponds to the inclusion of a restoring force into the 
equation of motion for the electrons within the material 
[22]. The dielectric function is then 



e(r,a;) = 1 - e(z)- 



w^ 



■ UJ^ 



(55) 
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coefficients of the surface [26], just as our formulae (41) 
and (42) do. Thus, for an end result that can depend 
only on the surface's reflection coefficients, one necessar- 
ily gets the same expressions (41) and (42), as before, 
for the magnetic moment shifts. Likewise, the Feyn- 
man propagator of the Huttncr-Barnctt field theory con- 
structed in Refs. [23, 27] depends only on the same reflec- 
tion coefficients, thus leading to the same conclusion that 
such a field theoretical approach necessarily reproduces 
Eqs. (41) and (42) for the magnetic moment shift. 

The continued validity of Eqs. (41) and (42) for a dis- 
persive dielectric is facilitated by the fact the integra- 
tion path C in the complex kz plane is not affected by 
the introduction of the transverse optical resonance ujt- 
The complex k^ plane for a dispersive dielectric model is 
shown in Fig. 6 with the abbreviations 
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FIG. 6. Complex k^ plane for dispersive dielectric surface. 



The dispersion relation for the surface polariton turns 
from Eq. (47) into 



1 



^sp = fc|| + ^H + 



l)-\ ki-ki 



'll'^T 



,2\2 



(56) 
The quantization of the electromagnetic field in terms of 
normal modes is now hindered by the fact that the field 
equation 



^z,± 
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*|| (jJrp ■ 
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While the structure of the complex k^ plane is consider- 
ably more complicated for the dispersive dielectric model 
as opposed to the plasma model, none of the additional 
cuts or poles interfere with the integration path C, and 
thus Eqs. (41) and (42) can be applied as they are. 



IV. EVALUATING THE MAGNETIC MOMENT 
SHIFT 



V X [V X AQ{r,u)]- 



e(z)a;2 

(57) 
cannot be written as a Hermitian eigenvalue problem, 
which would have guaranteed the orthogonality and com- 
pleteness of the modes. Thus, for a first-principles deriva- 
tion of the magnetic moment shift for this kind of surface 
one would need to include both dispersion and absorp- 
tion into the model and construct a Huttner-Barnctt- 
type field theory for the electromagnetic field interacting 
with the dielectric medium (for a suitable formulation see 
e.g. [23]). 

An alternative approach [24] is to use a Lifshitz-stylc 
method [25] and write the expectation values of squares 
of electromagnetic field operators, as in Eq. (21), in terms 
of a Green's tensor with an arbitrary permittivity e(r, lu). 
This tensor turns out to depend only on the refiection 



In order to calculate the magnetic moment shifts (41) 
and (42) for the various models of the surface, we di- 
vide up the integrals into the contribution from the two 
straight lines on cither side of the cut, kz = {—ioo, —iku] 



and kz = [— ifc 



II- 
and the contribution from the 

small circle around kz = — ifc|| . The integral around this 
small circle proves awkward to evaluate, so we subtract 
that point and consider it separately. This separate in- 
tegral turns out to be elementary; the result appears be- 
low as the terms outside the integrals. We write the full 
results in terms of the complex frequency ^ = iw and of 
?7 = kz/uj, which is the cosine of the complex angle of inci- 
dence. For isotropic media, this means that the dielectric 
function e is a function only of ^ and the material param- 
eters of the surface, like n, ujp, and wr- The magnetic 
moment shifts for the magnetic field Bq perpendicular 
and parallel to the interface are 
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(58) 
(59) 



with 
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TE 



^-V(g(e)-i)+^ 
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•TM 



v<o - v(e(e)-i) + ^2 



r;e(0 + V(KO - 1) + '/' 



(60) 



r 



Wc note that there is no need for the subtraction from 
i?TE of RTE{kz — ^ ~*^||) because that is zero. 

As shown later, for certain dielectric functions there 
may arise problems due to non-commutation between 
limits of physical parameters and the limit the k^ — t- 
— iA;||, which corresponds to the static limit w ^- 0, and 
this will prove to be crucial to the analysis and compar- 
ison of the results for the various models of the surface. 
The subtraction of the point at k^ — > —ikn is, in these 
variables, the subtraction of the point {^ — >■ 0,?7 — )■ oo}. 
Taking this two-dimensional limit risks the obvious pit- 
fall of potential non-commutation of the ^ and rj limits. 
This complication arises with the TE part of the inte- 
grals for the plasma surface. We emphasize that this is 
not an artefact of using these variables — even in an al- 
ternative formulation with {^ — > 0, r/ — > cx)} expressed as 
a single-variable limit, the issue manifests itself, though 
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in a different way. Fortunately, this one problematic case 
can be dealt with differently, as detailed in appendix C. 
Thus, the above formulae deliver the magnetic moment 
shift for a specified dielectric function e(r,a;), as long as 
the limits ^ — ?> and 77 — )- c» of the reflection coefficients 
commute. 

If one carries out the corresponding calculation in ci- 
ther a Huttner-Barnett field theoretical approach, as in 
Refs. [23, 27], or a noise-current approach [24], then the 
magnetic moment shift is more naturally expressed not 
as an integral over {k^, fcy), as in Eqs. (41) and (42), but 
as one over (a;,fc||). Hence, for ease of comparison with 
such approaches, we now proceed to derive alternative 
expressions for the magnetic moment shifts. Carefully 
considering the complex variable transformation from kz 
to uj one can transform Eqs. (41) and (42) to 
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e3 f°° f kn 

e3 f°° f ku 
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2u^ R 



»-TE 



fcjl + W^ ) rL 



„2ik. 



TM 
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(61) 
(62) 



r 



with the transformed path C" as shown in Fig. 7 and 
kz = —i^/k?. — oj-^ for |a;| > fc|| along C". These formulae 
can be simplified by subtracting the l/a;^ term of the 



integrand's Laurent expansion around a; = and treating 
it separately, when it integrates to zero along C". Then, 
using the fact that the integrands arc even in uj and re- 
writing the integrals in terms of the complex frequency 
^ = io;, we obtain 
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dfc|| / d^^l . 



Wn^m^J, *"io '^^W 
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(64) 
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FIG. 7. Integration path C" in the complex uj plane. 



and i?^]y[(0) = -Rtm(^ ~ 0). This result can also be 
obtained from a noise-eurrent approach [24]. 



The results in Eqs. (63) and (64) arc completely equiv- 
alent to those in Eqs. (58) and (59). Which ones are 
preferable depends on the particular model, but gener- 
ally speaking the difficulty in evaluating either is about 
the same for most dielectric functions e(^ = icu). 



V. RESULTS 



A. Non-Dispersive 



Inserting the dielectric function e(a;) = n? into 
Eqs. (58)-(60), one can evaluate the integrals exactly. 
The results for the magnetic moment shifts in a mag- 
netic field perpendicular and parallel to the surface are 
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(71- l)Vl + n' 
+ {n~ 1)71 



(66) 



If we expand these shifts in a series for large values of the 
refractive index 71 we obtain 
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0(1/71) 



(68) 

which shows that the perfect-rcfiector limit does not exist 
and the shift diverges for n —> 00. This is evidently un- 
physical, as it implies that the magnetic moment could be 
increased arbitrarily by increasing the refractive index 71 
of the surface. However, as we shall show below, this is in 
fact a misconception and it is that model of a dispersion- 
less medium which is unphysical. Detailed comparison 
with the results for the magnetic moment shifts near a 
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perfect reflector (cf. Eq. (7.12) of [9]) reveals the rather 
curious fact that the next-to-leading terms independent 
of 71 in Eqs. (67) and (68), if taken on their own, do in 
fact reproduce the results of the perfect-reflector case. It 
is easy to check that this is not a simple calculational 
error: if one takes the limit 71 —> 00 in the reflection coef- 
ficients (60) first and evaluates the equivalent of integrals 
(58) and (59) afterwards (as shown in Appendix B), one 
reproduces the results of Ref. [9]. The calculation also 
reveals the mathematical origin of the discrepancy: the 
TE reflection coefficient differs depending on whether ei- 
ther the perfect reflector limit n — > cx) or the static limit 
{^ — > 0,77 — > 00}, or kz — )■ — ifcii, is taken first. This 
indicates that the correct magnetic moment shift close 
to a specific material can only be obtained if the model 
chosen for the surface correctly reproduces the true low 
frequency behaviour of the dielectric susceptibility of the 
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material. Wc shall elaborate on this point later on. 



B. Plasma Surface 

Inserting the dielectrie function (43) into the reflection 
coefficients (60), we find that the 77 — >■ oo and ^ — >■ 
limits of i?TE do not commute. This means we cannot 
use the TE parts of Eqs. (58) and (59) for this particular 
model. However, for the TE parts of these integrals we 
can go back to the stage before we deformed the contour 
and carry out the integration along the path C, because 
the TE reflection coefficients are very simple so that the 

I 



integrals are unproblematic to calculate. As illustrated 
in Fig. 5, the contour C before deformation runs straight 
along the k^ axis, passing under the cut. We show in 
Appendix C how to calculate the integrals along the path 
C for the TE contributions. The contribution from the 
TM modes can be evaluated from Eqs. (58) and (59). To 
calculate the integrals we first replace the integration over 
77 by one over k, = ^-q. In the resulting two-dimensional 
integral we can then change the order of integration and 
go from Jg d^ jc dn to Jp '^'^/q d^- The integration 
over ^ is elementary, and we are left with just a one- 
dimensional integral over k. Simplifying the latter by 
scaling k to s = K/up we obtain 






167r2m3 1 ~ 4^^ '^P 



Wn^m^ 



;2+^P 



\g^2su:,. l+t\a) 



■ [2t{s) - (1 + 2t^(s)) arccot(t(s))] 



(69) 



"dse^-P- l + t\s) 



t2(s)[2~fi2(s)]3/2 

[3t(s) - (5 + 3f^(s)) arccot(i(s))] I, (70) 



i2(s)[2 + f2(5)]3/2 



with the abbreviation 



t{s) 




(71) 



Note that we have included the constant terms from 
Eqs. (58) and (59) since these originate from the TM 
reflection coefficient. The remaining integral over s is 
sufficiently complicated that it is best done numerically, 
which is easy since the exponential ensures fast conver- 
gence. We add to this the TE contribution calculated in 
Appendix C and show the full results in Fig. 8. 

The WpZ — > 00 asymptotics of Eqs. (69) and (70) 
is straightforward to work out by applying Watson's 
lemma. One finds that the s integrals give terms of order 
l/{ijjpz'^) and are thus negligible compared to the l/2;2 




FIG. 8. Comparison of plasma (solid line) and perfect re- 
flector results (dashed line). (Both axes are in dimensionless 
natural units.) 



I 

leading-order contributions in the first summands and 
from the TE part (worked out in App. C). The LOpZ — > 00 
limit of the overall result reproduces the perfect reflec- 
tor result (and the n independent terms in Eqs. (67) and 
(68)), as shown in App. C. This makes sense physically, as 
the lack of evanescent modes on the vacuum side means 
that for Wp — > 00 the plasma model becomes equivalent 
to the perfect reflector, thus reproducing its results. 

By contrast, the plasma model does not reproduce the 
results of the non-dispersive dielectric for large n. The 
consideration of a dispersive dielectric model in the fol- 
lowing section will shed further light onto this issue and 
point to the origin of this discrepancy in the very different 
low-frequency behaviour of the electromagnetic response 
of conductors and insulators. 

We also note that for the plasma model the Wp — > 
and 77 — > 00 limits of the TM reflection coefficient do 
not commute, causing the magnetic moment shift to di- 
verge as ojp — > 0, when it should clearly be zero at this 
point. This arises because the static limit of the dielec- 
tric constant is ill-defined for Wp — >■ and has already 
been discussed in detail in the context of the mass shift 
of an electron near a plasma surface [14]. Of course, if 
we take Wp — >■ before carrying out any of the integrals 
over photon wave vectors, there is no such problem and 
the shift vanishes as expected. 

We have shown above that the shift for the non- 
dispersive dielectric has a distance dependence of l/2;2 
for all z and all n, cf. Eqs. (65) and (66). However, for the 
plasma model we find that for small distances, i.e. small 
IwpZJ, the shift varies as l/z^. The TE part of the shift 
contributes only logarithmically at small IwpZJ, as shown 
in Eq. (Cll). So, the short-distance asymptotics of the 
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shift is dominated by the TM parts of Eqs. (58) and (59), 
which we have already simphfied into Eqs. (69) and (70). 
To find the small \u)pz\ asymptotics we scale the integra- 
tion variable s to a new variable equalling sLOpZ and then 
expand for small |wp2|. The resulting series may then be 
integrated term by term and gives 



Non-dispcrsive , ^ ' 



^^i^{\u:pz\ « 1) = ^-j-^ 



1 1 
16\/27r WpZ 

5 1 



+ 0{ujpz) 



(72) 



+ 0{LJpZ) 



(73) 



The leading 1/z^ dependence seen here arises because 
at small distances the interaction between the electron 
and the surface is dominated by electrostatic interaction 
of the electron with the surface plasmon. Therefore one 
should be able to derive the 1/z^ term by considering 
the surface plasmon part of the mode functions Eq. (49) 
alone. To this end, using the dielectric function (43) and 
the dispersion relation (48), we write the norming func- 
tion p(fc||), Eq. (50), and the imaginary z component of 
the wave vector k = ikz as functions of A:|| and ujp and 
then substitute them into Eqs. (53) and (54). Chang- 
ing variables such that the z dependence is taken out of 



the exponential and expanding for small 



we find 



the leading term for the perpendicular case given by the 
trivial integral. 



16to^v27i 
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(74) 

which is in agreement with the 1/z^ term in Eq. (72). 
For the parallel case the corresponding calculation repro- 
duces Eq. (73), as expected. Thus the 1/z^ dependence of 
the shift at short distances does indeed originate entirely 
from the surface plasmon part of the mode functions, as 
observed in Ref. [28] for the magnetic energy shift of a 
neutral atom interacting with a plasma surface. 

For large distances, the dependence of the magnetic 
moment shift remains l/^^, as shown in App. C. Since the 
inverse of the plasma frequency ujp is the only length scale 
in the plasma model large and small distance regimes are 
defined by \ujpz\ ^ 1 and \u!pz\ <^ 1. respectively, and 
hence to speak of large distances is the same as speaking 
of large ujp. 



C. Dispersive Dielectric 

A dispersive dielectric has features in common with the 
non-dispersive model, e.g. that modes originating from 
within the material can be evanescent on the vacuum 
side. Likewise the TE reflection coefficient exhibits the 
same kind of problem with non-commuting limits: the 
limit of a large dielectric response, which is now described 
by ojp — > oo, is not interchangeable with the static limit 




FIG. 9. Magnetic moment shift for dispersive and non- 
dispersive dielectric models as a function of static suscepti- 
bility, for the case of the magnetic field Bo perpendicular to 
the surface and \ujtz\ = {0.01,0.015,0.02.} 
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FIG. 10. Magnetic moment shift for dispersive dielectric and 
plasma models as a function of scaled distance \ujpz\ from the 
surface, for the case of the magnetic field Bo perpendicular 
to the surface and LJr/i^p = {0.02, 0.04, 0.06.} 



a; — 7> 0. Thus for large e we expect the dispersive dielec- 
tric to give similar results to the non-dispersive dielectric. 
We evaluate Eqs. (58) and (59) numerically, and find a 
peak in the magnetic moment shift relative to the perfect- 
reflector result. To facilitate the discussion of this peak 
and the comparison of different models, we now choose 
to write the dielectric function in terms of the static limit 
of the dielectric susceptibility. 



X(0) = e(0) - 1 - c.2/c,2 



We find peaks in Afi± and A/X|| at \/x(0) « 2, with 
the height of the peak being inversely proportional to 
lotz, as shown in Fig. 9 for the case where the external 
magnetic field Bq is perpendicular to the interface. We 
also plot the corresponding shift for the non-dispersive 
case, where x(0)nondisp = n^ — 1. If we were to con- 
tinue the plot to very large values of x(0), the graphs 
for the two models would very slowly converge into one 
linearly-rising curve. By contrast, the result for the per- 
fect reflector, also shown in Fig. 9, is much smaller and 
has the opposite sign. 

The peak appears if the choice of parameters is such 
that \u!tz\ < 0.07 for Bq perpendicular, and < 0.25 for 
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Bo parallel to the surface. For smaller values of Iwt^I , the 
peak moves closer to \/x(0) = uip/uiT ~ 2, and increases 
in height. To gauge the enhancement that dispersion 
brings to the shift we calculate the ratio of the height 
of the dispersive peak to the non-dispersive result at the 
same x(0), and find 



A/i 



_Ldisp 



^M-Lnondisp 



30.3 eVnm 

\u)tz\ 



^A^lldisp 



Ami 



nondisp 



81.6 eVnm 

\ujtz\ 



(75) 

A typical value for the frequency uit in a metal is on the 
order of a few eV (see, for example, [29]), meaning that 
a significant enhancement relative to the non-dispersive 
case would be observed only at extremely small distances 
z. However, restricting oneself to considering the prop- 
erties of only elemental solids would be short-sighted. 
Structures engineered on the nanoscale can have trans- 
verse resonance frequencies lot significantly smaller than 
any ordinary material — examples include an InSb semi- 
conductor grating with lot (and LOp) in the range of a few 
mcV [30] . These types of materials are at a focal point of 
strong contemporary interest in low-frequency plasmon- 
ics. With appropriate assumptions about the approxi- 
mation of a part of such a structure as a planar surface, 
we find that for distances z of a few tens of nanometres 
one may get an enhancement factor on the order of 10'^ 
relative to the non-dispersive case. 

The apparent problem of the behaviour of the non- 
dispersive result in the limit of large refractive index, 
n — >■ 00, can be clarified by comparing it with the be- 
haviour of the dispersive shift at large x(0)- In this 
regime the shift for the dispersive dielectric model be- 
comes linear in ^/xW) aud agrees with the non-dispersive 
results; so for large x(0) the two models are equiva- 
lent. The crucial additional observation is to note that 
for a non-dispersive dielectric with large x(0) we have 
x(0) ~ n^, which is to say that a large the refractive index 
necessarily implies a large static susceptibility. There- 
fore, in the non-dispersive model one cannot sensibly 
make a distinction between an arbitrarily large refrac- 
tive index and an arbitrarily large static susceptibility. 
Investigation of the dispersive dielectric has shown that 
the latter interpretation is the correct one — the mag- 
netic moment shift grows with increasing static suscepti- 
bility, but an arbitrarily large static susceptibility is, of 
course, physically impossible. So while the shift in the 
non-dispersive case does indeed increase without bound 
as the refractive index n is increased, this is not due any 
problem with the calculation, but is in fact the result 
of the static susceptibility growing without bound and 
an inevitable consequence of the unrealistic exclusion of 
dispersion from the model. 

Consideration of the shifts in terms of the static sus- 
ceptibility also emphasizes the close relationship between 
plasma and perfect reflector models. In both of these 
models the static susceptibility is infinite right from the 
start, which means that their results do agree in the limit 

Wo — >■ OO. 
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FIG. 11. Commutation properties of the various models dis- 
cussed in the text. Each arrow indicates a limit which takes 
one dielectric function to another. Solid (dashed) arrows in- 
dicate a limit that, when applied the reflection coefficients, 
commutes (does not commute) with the limit kz -^ ~*^|| • 
The consequence of this is that the magnetic moment results 
for two models connected by solid (dashed) arrows are (are 
not) obtainable as limiting cases of one another. 



The differences between the four models discussed 
above very clearly show that in order to predict the mag- 
netic moment shift for a given set-up, one must choose 
a model which is physically appropriate for the low- 
frequency behaviour of electromagnetic response of the 
material at hand. In other words, it matters whether 
the material is a conductor or an insulator. These two 
classes of material are not obtainable as limiting cases 
of each other because the conductor models ignore the 
existence of evanescent modes (which is a direct conse- 
quence of their static susceptibilities being infinite) . The 
calculations for each class of model diverge from each 
other because of non-commutation of a variety of limits 
of the reflection coefficient, namely between the static 
limit (fcz — > —ik\\) and whichever limit we have to take 
in order to compare models. For example, we note that 
the 71 — > oo and k^ — > ^^^jl hmits of the non-dispersive 
TE reflection coefficient do not commute, which leads to 
the n — >■ oo limit of the result for a non-dispersive di- 
electric to disagree with the perfect refiector result. A 
further important example is that the limit of vanishing 
transverse resonance frequency, wt — > 0, and the static 
limit kz — > — i^ii of the dispersive TE refiection coeffi- 
cient do not commute, which means that taking ujt — >■ 
(x(0) ~^ oo) in the dispersive dielectric results will not 
reproduce the plasma results, while naive comparison of 
the dielectric functions (43) and (55) suggests that they 
should. Since this plays such an important roles for the 
physical interpretation of the results, we summarize the 
commutation (or lack thereof) between the various limits 
of the refiection coefficients in Fig. 11. 

For the plasma model, we found a Xj z^ dependence of 
the magnetic moment shift at small distances, i.e. small 



17 



Iwpzj, and that the leading 1/z^ term can be found either 
by determining the asymptotics of the complete shift, or 
by considering only the part due to the interaction with 
just surface plasmons. The asymptotics of the integrals 
for the shift in the dispersive dielectric case are too awk- 
ward to analyse directly. Instead we give the results one 
obtains by considering only the interaction with the sur- 
face polariton, i.e. by using Eqs. (55) and (56) in Eqs. (53) 
and (54): 



A^l^{\uJpz\ < 1) 



A/X||(|a;pz| < 1) 



64:Tr\/2in'^z^ 



be-' 



2w|,- 



1287rV2m3z3 



2w|,- 



(76) 



(77) 



Surprisingly, for these short-distances tUpZ <C 1, we find 
that the ujt -^ limits of Eqs. (76) and (77) do agree 
with the corresponding results for the plasma, Eqs. (72) 
and (73), unlike the results for general distances ujpZ > 1. 
This is because these results depends only on the surface 
plasmon part of the mode functions and electrostatic in- 
teractions, but there is no refiection of travelling photon 
modes, and hence any non-commutation of limits in the 
reflection coefficient does not come into play. 

At large distances ojpZ :$> 1 the shift decreases as l/z^, 
which is obvious when applying Watson's lemma to the 
integrals in Eqs. (58) and (59) or, more conveniently, to 
those in Eqs. (63) and (64). In fact, its asymptotics is 
given by the same expressions as in Eqs. (65) and (66) 
except with n replaced by y'l -I- {uJp/ujtY, i.e. by the 
square root of the dielectric function at zero frequency, 
which is expected because at large distances only the 
static electromagnetic response of the surface matters 
and dispersion plays no role to leading order. Figure 
10 shows the shift for dispersive dielectric and plasma 
models as a function of distance. 



VI. EXPERIMENTAL RELEVANCE 

Expressing magnetic moment shifts as relative shifts 
A^/^ to the Dirac magnetic moment jjl ~ e/277i, we have 
for the perpendicular component of the non-dispersive 
shift in SI units; 



^A^±nondisp 



c'eo \&T^ m^z- 



r/W 



10""nm2 



(78) 



where f{n) is the remaining part of Eq. (65), and is of 
order unity. For a distance z w Inm, Eqs. (75) (and the 
discussions following them) show that the enhancement 
due to the inclusion of dispersion can be of order 10"* 
under favourable conditions. Thus, we have a magnetic 
moment shift of up to one part in 10^. The current ex- 
perimental accuracy for 17/2 in free space is on the order 
of one part in 10^^ [31], so that the shift calculated here 



would compare very favourably to this. As the distance 
increases to the order of a micron the effect decreases to- 
wards the limits of current experimental accuracy. For 
example, an electron 0.1//m away from the same surface 
as above would have its magnetic moment shifted by only 
one part in 10^^. 

This leads one to ask if the current best techniques for 
measuring the g factor would be suitable for making a 
measurement of the surface dependent shift of the mag- 
netic moment. Since one of the sticking points in such ex- 
periments is that accurate measurement of the externally 
applied magnetic field Bq is mostly impossible, g-factor 
experiments usually do not measure the magnetic mo- 
ment directly, but instead they find its ratio to either a 
known magnetic moment, or to the cyclotron frequency 
of the particle under consideration. In case of the lat- 
ter for surface-dependent magnetic moments shifts one 
would need to take into account the shift in cyclotron fre- 
quency of a particle near a surface, which arises due to the 
position-dependent self-energy of the particle [9, 13, 14]. 
Crucially, the leading term of the surface-dependent cy- 
clotron frequency shift is of order a/{mz) and thus much 
bigger than the magnetic moment shift which is of order 
a/(mzY (cf. Eq. (78)). So, an experiment which adapts 
the techniques used for measuring the free-space g fac- 
tor to find its surface dependant part would effectively be 
measuring the change due to the surface in its self-energy, 
not in its magnetic moment. While direct experimental 
confirmation of a shift in the self-energy would, of course, 
be interesting in its own right, its existence represents a 
significant obstacle to isolation and observation of the 
magnetic moment shift. 



VII. SUMMARY AND CONCLUSIONS 

We have shown that the magnetic moment shift near a 
non-dispersive imperfectly reflecting surface has notable 
differences from the corresponding shift near a perfectly 
reflecting surface, differing by orders of magnitude and 
in some cases even in its sign. The inclusion of dispersion 
can significantly modify the magnitude of the effect, and 
that this modification can be tuned by choice of mate- 
rial. We have given a general formula for the calculation 
of magnetic moment shifts, into which one can simply 
insert the relevant dielectric function and evaluate the 
integrals to obtain the shift as dependent on the distance 
from the surface and various parameters characterizing 
the electromagnetic response of the surface. 

The work presented here has also considerably ex- 
tended the horizon of the traditional theory of cavity 
quantum electrodynamics. While the effect that causes 
the shift of the magnetic moment of an electron close to 
a surface is essentially the same as the one that causes 
the Casimir-Polder shift of energy levels in neutral atoms 
close to a surface, the theory required for the magnetic 
moment shift is considerably more involved. On the one 
hand, this is because the nature of the spin requires a 
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rclativistic treatment, but on the other hand, the much 
more subtle but ultimately more important reason is 
that dynamics of the system is qualitatively more compli- 
cated: a neutral atom interacting with a surface involves 
the exchange of virtual excitations between just the pho- 
ton field and the atomic electron, whereas the shift of 
magnetic moment of an electron near a surface involves 
virtual excitations being exchanged in a triangle between 
the photon field, the centre-of-mass motion of the elec- 
tron, and its spin states. 



In terms of these the canonical momenta are then given 
by 



^y=Py 2^^ = MbR + bj^) 



i.-\ 



so that the Hamiltonian reads 
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Appendix A: Schrodinger and Dirac equations for a 
particle in a constant magnetic field 



In order to derive the eigenstates of the Dirac equation 
for a particle in a constant classical magnetic field Bq one 
first solves the corresponding Schrodinger problem and 
then uses its solutions to generate the Dirac eigenstates 
[17]. The Schrodinger Hamiltonian for a charged particle 
moving in a constant magnetic field Bq = Bqz is 



Hs 



iP. + ^y)' , iPy-^x? , pI 



2m 



2m 



2m 



(Al) 



One can reduce this to a harmonic oscillator by introduc- 
ing annihilation and creation operators and rewriting the 
positions and momenta in terms of those. 



i=^(6, + 5t) p,, = !^(6t_5j (A2) 

y = j^^K+H) py-'-^iK-K) (A3) 



This is also the reason for why there is no simple reso- 
nance effect when the spin flip frequency in a magnetic 
field coincides with an absorptio n resonan ce in the mate- 
rial of the surface, where (3q = y^—eBo/2. (Note that we 

use e = — |e|). The operators b.^ and 6j, are then combined 
to form creation and annihilation operators for right and 
left-circular quanta 

^fl = —^ik- i'by) , 6^ = -^{bl+ ibl) , 
K = ^{K+^by), 6[ = -L(6t_,St). 
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■p 
2m 
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Thus the Hamiltonian is equivalent to a harmonic oscilla- 
tor of right-circular excitations and possesses infinite de- 
generacy with respect to the left-circular quanta. Eigen- 
states |i^) of the Schrodinger Hamiltonian Hs can there- 
fore be generated by repeated application of the creation 
operator 6^ to the ground state \iy = 0) which is defined 
by b^\i^ = 0) = 0. 

We can now use the Schrodinger eigenstates to derive 
the corresponding Dirac eigenstates. Following [17], we 
start by noting that eigenfunctions of the Dirac equation 

(a • 7r + fim)ijj = Hqi/j = E,,^ (A5) 

may be obtained from solutions of 

{Hi - El)X = {Ho - E,)iHo + E,)X = . (A6) 

Evidently, if a state X satisfies the above equation, then 

i' = (Ho + E,)X (A7) 

is a solution of the Dirac eigenvalue problem, Eq. (A5). 
To find the eigenvalues E'^ of H^ wc calculate iJg , which 
on account of(Q:-7r)^=7r^ — ecrzBo and the anticom- 
mutation properties of the Dirac matrices, reads 



Hq = TV — CCTzEq + 771 



(A8) 



Therefore we can relate Hq to the Schrodinger Hamilto- 
nian Hs^ 



Hq = 2mHs — ecTz-Bo + m, 



(A9) 



and then the eigenvalues E^ of H^ follow from Eq. (A4) 
and the eigenvalues s of the spin operator S^ = crz/2, 



Et =m^+pi-2eBo U + s 



(AlO) 



We now choose the states X in such a way that they dis- 
tinguish spin-up and spin-down states, and particle and 
anti-particle states, i.e. we choose them to be eigenfunc- 
tions of az with eigenvalues s = ±1/2, and of /3 = 70 
with eigenvalues 1 for a particle and —1 for an antipar- 
ticle. Equation (A9) implies that the Dirac eigenstates 
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can be expressed in terms of a product state of the non- 
relativistic eigenstates \v) and the spin state |s), which 
we choose to write as \v,s) = \i') ®\s), 



I*. 



if + El, 



^2E,AE, + m) 



|i^) x^'^'^^ for s = ±1/2, (All) 



where x*'^^^ = (1,0,0,0), x^^^'^ = (0,1,0,0). For an- 
tiparticle eigenstates the negative root of (AlO) ap- 
phes, the normahzation factor in the denominator of 
Eq. (All) turns into yJ—2E,j{—Ey + m), and we use 

x(^)^ = (o,o,i,o),x(^^^t = (o,o,o,i). 

For calculations it is expedient to express momentum 
components in terms of 

7r+ = TTj; + illy = 2i/3o6|j. , 
7r_ = TTj^ — iiTy = —2ij3Qbf{ 

Thus, for a general vector Q we have 



Q • TT = z/^oQ-^R - i/3oQ+bR + Q^p^ 



with 



Q+ = Qx + iQy 
^/_ ^x ^^y 

In the main text we also need the matrix elements of 
the position operator. Noting that [32] 

^ = ^i^R + bL+b^R + bl) ^xo + --{bR + b^j,) 

V = 7^{bR -bL-Vj, + bl) =2/0 + ^{bn - 6^) , 
ZPQ zpo 

we have 
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1 

^Po 



(A15) 



not apply to the perfect reflector which has a TE reflec- 
tion coefficient of —1 at all frequencies, so that Eqs. (58) 
and (59) cannot be used for calculating the magnetic mo- 
ment shift near a perfect reflector. Instead, we go one 
step back and start from Eqs. (41) and (42). At flrst 
glance it appears that the branch cut due to /cf might 
meet that due to fc when the limit n — >■ 00 is taken. How- 
ever, in fact the branch cut due to /cf disappears since 
the perfect reflector excludes all right-incident modes, so 
that that /cf does not appear in the integrand. Therefore 
the contour of integration simply runs straight along the 
kz axis. This can also be seen formally by using the inte- 
gral over C", taking the limit n — >■ 00 (which eliminates 
the branch cut) and then deforming the contour back up 
to the real axis. Either way one obtains 
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where we have carried out the k^ integration flrst and 
used 



/fell e2*'=-^ 



dk^ 
- (A:2 + k^ 



3/2 
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Eqs. (Bl) and (B2) reproduce the shifts calculated by [9], 
and the n- independent terms in Eqs. (67) and (68). This 
calculation shows another way of looking at the funda- 
mental disparity between the perfect reflector and non- 
dispersive dielectric models: the branch cut due to fcf 
meets that due to k if n is taken to be flnitc at the start 
of the calculation and then made inflnite at the end, while 
if n is infinite from the start the branch cut never appears 
in the flrst place. 



Appendix B: Reproduction of perfect mirror result 

In the derivation of Eqs. (58)and (59) it was argued 
that the subtraction and re-addition of the point {^ — )• 
0,?7 — >■ 00} from the TE reflection coefficients would not 
be necessary since R^^ is zero at this point for any phys- 
ically reasonable dielectric function. However, this does 



Appendix C: TE Part of Plasma Surface 

Using Eq. (41) but with the integration contour C" 
changed back to C for the plasma model before deforma- 
tion (cf. Fig. 5), one can write the TE part of the shift 
for the plasma surface as 



A/x±, 



TE 



32T:^(jj'im^ 



dk\\ 



dk- 



hj\\\ ZifZu rZ y, j 

(fcjj + fc2)3/2 



2ki 



2k, 



fc? - W2 ) g2zfc. 
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(CI) 



where the reflection coeflicient has been written out exphcitly using Eq. (45) and the branch cut is taken to be 
between the branch points at ±Wp. Care must be taken when evaluating this integral due to the physical requirement 
sgn{kz) = sgn(fcf) from refraction; so we outline the calculation here. First, we note that the order of integration 
matters, and that the integral is convergent only if the k^ integration is carried out first. To circumvent this problem 
we introduce a cutoff A on the fcii integral. This improves the convergence of the double integral so that we are allowed 
to interchange the order of integrations. The fc" integral can then be calculated exactly and gives: 



Am±, 



TE 






dk, (2A - 5|A;,|) (2kl - co^ ~ 2k,,^kl-ul) e^.fc,. ^ ^^^/^^ 



(C2) 



We first consider the term with the square root. For the contribution from the region \kz\ > tOp-, we have (omitting 
the overall constants): 



dkz 



dk. 



„2ik. 



= (2A-5|fc,|)(-2fc,^fc2_^2^ 



(C3) 



Noting that k,^ k1 — iJi is even in k, because of the which means that the factor k. 



iJi — k\ is now odd in 



physical constraint sgn(fcz) = sgn^/fc^ — w^, we can sim- 
plify this to: 



/■OO 

2 / dk. 



,cos(2fc,z)(2A-5|A:,|) (-2fc,y^fc2 „ ^2^ 



Next we consider the region jfe^j < Wp, where fcf 



fc2 — Wp is imaginary. As illustrated in Fig. 5, the inte- 
gration path runs along one and the same side of the cut, 
I 



kz- Applying the constraint sgn(fc2) = sgn^/fc^ — uSi, to 

the vicinity of k, ~ Wp, we are directed to choosing the 
complex sheet of the square root such that in the lower 

half-plane 



A: 2 — a; 2 = — i 



tions lead us to write the integral analogous to Eq. (C4) 
but from the region | /c^ | < Wp as 



(C4) 2 / " dk, (2A - S/c^) {-2k, Jul - kfj sin(2fc^z) (C5) 



The rest of Eq. (C2) is a trivial integral, and combining 
this with Eqs. (C4) and (C5) gives: 



gJ I r^p 

.f. , 3 lim S — / dk. 



(2/fc2 _ a;2) cos.{2k,z) - 2k, Jul - B, ii\n{2k,z) (2A - bk,) 
+ ^ f dk, cos{2k,z) (2fc2 - a;2 _ 2k,^k^, - ulj (2A - 5k,) i . (C6) 

p J LJp J 

The integrals proportional to A give expressions with the Bessel function J2{2ujpz), sin(2a;pz), and cos(2a;pz), but all 
together they conspire to add up to zero. Defining 



Ite = ^7-; / dk,k. 



'^p I JO 



we therefore have 



(2fc2-w2) cos{2k,z)-2k,Jujl - B, sin(2fc^; 



r 



dk, k, cos{2k,z) ((2fc2 - cj2) - 2k,^kl-Lol^ 

(C7) 



^p"^ pO 

Anj_^TE = -TT^ — ^T-TE and I^[i\\,te = ir^ — -T-te 



where the case for Bq parallel to the interface has been 
evaluated in exactly the same way. The integral T-te may 



IGir'^m^ 



2rir,-]0 



STT^m 



(C8) 
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be evaluated analytically; one finds 



X^ 



TE 



1 

4^ 



4za;^ ■nujpYi{~2ujpz) 



4z4o 



15 



2z 



'i'KY2{—2ujpz) TrH2{2ujpz) TTUjpH3{2ujpz) 



4^2 



4z2 



2z 



(C9) 



where y„ is the nth Bessel function of the second kind, 
and Hn is the nth Struvc function. This result displays 
the expected behaviour that lim^_>._oo Ite = 0, i.e. that 
there is no magnetic moment shift due to a surface that 
is infinitely far away. 

The 'perfect-mirror' limit of this object is 



lim Ite 



1 

4^ 



(CIO) 



which means that for the plasma surface the TE modes do 
not result in unlimited growth of the magnetic moment 



shift for large values of the permittivity, in contrast to 
what was observed for the non-dispersive dielectric in 
Eqs. (67) and (68). 

Inspection of Eqs. (69) and (70) in the limit WpZ — > oo 
shows that the s integrals give rise to terms proportional 
to 1/{ll>pZ^), whence only the l/z^ terms in front of them 
contribute to leading order. We then have in total 



A^±{llJp — > oo) 
A/i|[ {ujp — !> cxj) 



3 

I6n^m^ \ 4z2 4^2 

y2 



e 



e 



167r2m'^ \ 2z2 



32TT^m^z'^ 
327T^m'^z^ 



in agreement with the perfect- mirror limit, and also with 
the n- independent terms from Eqs. (67) and (68). 

We note that the asymptotics of Ite for small \u!pz\ 
are 



Ite{\c^pz\ « 1) = -yf [1 + 47 + 41n(-c^pz)] (Cll) 
where 7 is the Eulcr constant sa 0.577. 
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